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Abstract 

We consider Brownian particles with the ability to take up energy from the environment, 
to store it in an internal depot, and to convert internal energy into kinetic energy of motion. 
Provided a supercritical supply of energy, these particles are able to move in a "high velocity" 
or active mode, which allows them to move also against the gradient of an external potential. 
We investigate the critical energetic conditions of this self-driven motion for the case of a linear 
potential and a ratchet potential. In the latter case, we are able to find two different critical 
conversion rates for the internal energy, which describe the onset of a directed net current into 
the two different directions. The results of computer simulations are confirmed by analytical 
expressions for the critical parameters and the average velocity of the net current. Further, 
we investigate the influence of the asymmetry of the ratchet potential on the net current and 
estimate a critical value for the asymmetry in order to obtain a positive or negative net current. 

PACS numbers: 

05.40.Jc Brownian motion 

05. 45. -a Nonlinear dynamics and nonlinear dynamical systems 

05.60.-k Transport processes 

82.20.Mj Nonequilibrium kinetics 

87.10.+e General theory and mathematical aspects 
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1 Introduction 

The motion of a "simple" Brownian particle is due to fluctuations of the surrounding medium, i.e. 
the result of random impacts of the molecules or atoms of the liquid or gas, the particle is immersed 
in. This type of motion would be rather considered as passive motion, simply because the Brownian 
particle does not play an active part in this motion. On the other hand, already in physico-chemical 
systems a self-driven motion of particles can be found [1]. For instance, small solid particles floating 
on a liquid surface may produce a chemical substance at a spatially inhomogeneous rate, which 
locally changes the surface tension. This way the substance induces a net capillary force acting on 
the particle, which results in the particle's motion. 

Recent investigations on interacting self-driven particles show a broad variety of interesting phe- 
nomena, such as phase transitions and the emergence of self-ordered motion [2-4]. Here, the focus 
is on collective effects rather than on the origin of the particle's velocity; i.e. it is usually postulated 
that the particles move with a certain non-zero velocity. 

In this paper, we focus on the energetic aspects of self-driven motion in order to derive conditions for 
an active mode of motion. Active motion, as the name suggests, occurs under energy consumption, 
it is related to processes of energy storage and conversion into kinetic energy. For instance on the 
biological level, cells or simple microorganisms are capable of active, self-driven motion, which in 
many cases has been successfully described by stochastic equations [2,5-9]. 

In order to describe both the random aspects and the energetic aspects of active motion, we have 
introduced a model of active Brownian particles [10-15]. These are Brownian particles with the 
ability to take up energy from the environment, to store it in an internal depot and to convert 
internal energy to perform different activities, such as metabolism, motion, change of the environ- 
ment, or signal-response behavior. Since the focus in this paper is on the energetic aspects of active 
Brownian particles in a specific potential, possible changes of the environment are neglected here. 

A very simple mechanism to take-up the additional energy required for active motion, is the pump- 
ing of energy by space- dependent friction [11]. In this case, the friction coefficient 70 becomes a 
space-dependent function, j(r), which in a certain spatial range can be also negative. Inside this 
area the Brownian particle, instead of loosing energy because of dissipative processes, is pumped 
with energy, which in turn increases its velocity. Provided a supercritical supply of energy, the 
particle should be able to move with a velocity larger than the average thermal velocity. 

While such an approach will be able to model the spatially inhomogeneous supply of energy, it 
has the drawback not to consider processes of storage and conversion of energy. In fact, with only 
a space-dependent friction, the Brownian particle is instantaneously accelerated or slowed down, 
whereas e.g. biological systems have the capability to stretch their supply of energy over a certain 
time interval. 
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In order to develop a more realistic model of active motion, we have considered an internal energy 
depot for the Brownian particles [14, 15], which allows to store the taken-up energy in the internal 
depot, from where it can be converted e.g. into kinetic energy, namely for the acceleration of motion. 
Additionally, the internal dissipation of energy, due to storage and conversion (or metabolism in a 
biological context) can be considered. 

With these extensions, the Brownian particle becomes in fact a Brownian motor [16-18], which is 
fueled somewhere and then uses the stored energy with a certain efficiency [15] to move forward, 
also against external forces. Provided a supercritical supply of energy, we have found that the 
motion of active Brownian particles in the two-dimensional space can become rather complex [14]. 

In this paper we investigate the one-dimensional motion of an ensemble of Brownian particles with 
internal energy depot in piecewise linear potentials. One particular example is the ratchet potential, 
i.e. a periodic potential which lacks the reflection symmetry. Ratchet systems recently attracted 
much interest with respect to transport phenomena on the microscale, since they provide a mech- 
anism to transfer the undirected motion of Brownian particles into a directed motion. Hence, the 
term Brownian rectifiers [19] has been established. In order to reveal the microscopic mechanisms 
resulting in directed movement, different physical ratchet models have been proposed [20], such as 
forced thermal ratchets [21], or stochastic ratchets [22,23], or fluctuating ratchets [24,25]. 

The model discussed in this paper, aims to add a new perspective to this problem which is based 
on the idea of internal energy storage. In Sect. 2, we discuss the basic features of our model without 
specifying an external potential, and derive the equations of motion both for the general case and 
the overdamped limit. In particular, we point to the possible existence of a "high velocity" or active 
mode of motion for the Brownian particle, in addition to the usual passive mode of motion. In Sect. 
3, the necessary conditions for such an active mode of motion are investigated for the case of a 
linear potential. For the overdamped limit, we derive critical parameters of energy conversion, which 
allows the particle to move also against the direction of the external force, i.e. to move "uphill" 
the potential gradient. The results are applied to the motion of Brownian particles in a ratchet 
potential in Sect. 4. After investigating the deterministic motion of a single particle in Sect. 4.1, we 
discuss the energetic conditions for the establishment of a net current for the determinsitic motion 
of an ensemble of particles in Sect. 4.2. We show that, dependent on the conversion of internal into 
kinetic energy, a net current into different directions can be established. The critical parameters 
for the energy conversion and the resulting velocity of the net current, which are found by means 
of computer simulations, are compared with the analytical results obtained in Sect. 2. We also 
investigate the influence of the asymmetry of the ratchet potential on the establishment of the net 
current. 
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2 Equations of Motion of Pumped Brownian Dynamics 

The motion of simple Brownian particles in a space-dependent potential, U (r) can be described by 
the Langevin equation: 

r = v ; mv = —jo v — Vi7 (r) + T{t) (1) 

where 70 is the friction coefficient of the particle at position r, moving with velocity v. J-(t) is a 
stochastic force with strength D and a (5-correlated time dependence 

(Ht))=0; (F{t)F{t>))=2D5{t-t') (2) 

Using the fluctuation-dissipation theorem, we assume that the loss of energy resulting from friction, 
and the gain of energy resulting from the stochastic force, are compensated in the average, and D 
can be expressed as: 

D = k B T l0 (3) 
where T is the temperature and ks is the Boltzmann constant. 

In addition to the dynamics described above, the Brownian particles considered here are pumped 
with energy from the environment, which they can store in an internal depot. Further, internal 
energy can be converted into kinetic energy Considering also internal dissipation, the resulting 
balance equation for the internal energy depot, e, of a pumped Brownian particle is given by: 

j t e(t) = q(r) - c e(t) - d(v) e(t) (4) 

q(r) is the space-dependent pump rate of energy and c describes the internal dissipation assumed 
to be proportional to the depot energy. d(v) is the rate of conversion of internal into kinetic energy 
which should be a function of the actual velocity of the particle. A simple ansatz for d(v) reads: 

d(v) = d 2 v 2 (5) 

where d 2 > is the conversion rate of internal into kinetic energy. The energy conversion results in 
an additional acceleration of the Brownian particle in the direction of movement, expressed by the 
vector e v = v/v. Hence, the equation of motion for the pumped Brownian particles has to consider 
an additional driving force, di,e{t)v. In [14,15], we have postulated a stochastic equation, which is 
consistent with the Langevin Eq. (1): 

r = v ; mv + jqv + VU (r) = d2e(t)v + T{t) (6) 
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If we restrict ourselves to the one-dimensional case, i.e. the space coordinate is given by x and 
further assume m = 1 for the mass and q(r) = qo = const, for influx of energy, the dynamics for 
the pumped Brownian motion, eqs. (6), (4) can be specified as follows: 



x = v 



dU(x) 



v = -( 1 o-d 2 e(t))v-—^ + V2DZ(t) (7) 
e = qo — ce — d2V 2 e 
where £ is a stochastic force with white-noise fluctuations: 

(mm) = s(t-t') (8) 

In order to find a solution for the coupled Eqs. (7), let us now consider a relaxation of the dynamics 
on different time scales. If we assume that the velocity v(t) is changing much faster than the space 
coordinate x(t), for v(t), Eq. (7), a formal solution can be given: 



v{t) = u(0)exp j- 7o t + d 2 ^ e{t')dt'^ 
+ exp |- 7o t + d 2 j e{t')dt' 



jf exp j 7o t / - d 2 jf* e{t")dt" 



■VU + ^/2k B T 70 £(t')\ dt' (9) 

This solution however depends on the integrals over e(t), reflecting the influence of the energy depot 
on the velocity. If we further assume a fast relaxation of the depot, e(t), compared to the relaxation 
of the velocity v(t), the corresponding equation of Eq. (7) can be solved, and we obtain with the 
initial condition e(0) = the following quasistationary value for the energy depot: 

eo = , , 2 (10) 
c + d 2 v z 

which yields a possible maximum value of e™ ax = qo/c. 

The overdamped limit is obtained by considering a fast relaxation of the velocities, in which case 
the set of Eq. (7) can be further reduced to: 

.(*) = — i-7r + ^F«M (id 

7o - d 2 e Q dx 70 - d 2 e 

We note that, due to the dependence of eo on v 2 = x 2 , Eq. (11) is coupled to Eq. (10). Thus, the 
overdamped Eq. (11) could be also written in the form: 

■■<> " rf 2 J° . 2 ) x = -^ + y^fTo £(t) (12) 
c + d 2 x A J ox 
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Eq. (12) indicates a cubic equation for the velocities in the overdamped limit, i.e. the possible 
existence of non-trivial solutions for the stationary velocity. For the further discussion, we neglect 
the stochastic term in Eq. (12) and denote the stationary values of v(t) by vq(x). Further, the force 
resulting from the gradient of the potential, F(x) = — VJ7, is introduced. Then, Eq. (12) can be 
rewritten as: 



Depending on the value of F and in particular on the sign of the term ((70^2 — erf), Eq. (13) has either 
one or three real solutions for the stationary velocity, vq. The always existing solution expresses a 
direct response to the force in the form: 



This solution results from the analytic continuation of Stokes' law, vq = -F/70, which is valid for 
d.2 = 0. We will denote this solution as the "normal response" mode of motion, since the velocity 
v has the same direction as the force F resulting from the external potential U(x). 

As long as the supply of the energy depot is small, we will also name the normal mode as the passive 
mode, because the particle is simply driven by the external force. More interesting is the case of 
three stationary velocities, vq, which significantly depends on the (supercritical) influence of the 
energy depot. In this case which will be discussed in detail in the following section, the particle will 
be able to move in a "high velocity" or active mode of motion. For the one-dimensional motion, in 
the active mode only two different directions are possible, i.e. a motion into or against the direction 
of the force F. 

But already in the two-dimensional case there are infinitely different possibilities. This conclusion 
is of importance when discussing stochastic influences as will be done in a forthcoming paper [26]. 
In the one-dimensional case, the influence of noise is rather weak, because of the limited number of 
possible directions, but in the two- and three-dimensional case, the active motion of the particles is 
very sensitive to stochastic influences, which may determine the direction of motion in the active 
mode. 

3 Deterministic Motion in a Linear Potential 

3.1 Stationary Solutions 

In the following, we restrict the discussion to the one-dimensional, deterministic motion of the 
particle, corresponding to D = in Eq. (7). Further, we may assume that the force resulting from 
the gradient of the potential is constant. Then, we have the two coupled equations for v(t) and 



^270 Vq - d 2 Fv - (q d 2 



C7o) Vq = cF. 



(13) 



v (x) ~ F(x) 



(14) 
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e(t): 

v = -(~f -d 2 e(t)^v + F (15) 
e = qo — ce — d2V 2 e 

The stationary solutions of Eq. (15) are obtained from v = and e = 0: 

v = — ; e = 2 (16) 

7o - a 2 e c + d 2 v 5 

which lead to the known cubic polynom for the amount of the constant velocity i>o, Eq. (13): 

d 2 -f vl - d 2 Fvl - (q Q d 2 - C70H - cF = 0. (17) 
Here Vq +1 is defined as a vector |t>o| n t>o- 

For a first insight into the possible solutions of Eq. (17), we consider the force free motion, F = 0. 
Then Eq. (17) has either one or three stationary solutions which read: 



^-S^-k (18> 

The corresponding bifurcation diagram is shown in Fig. la which displays a typical fork bifurcation. 
The bifurcation point is given by: 

df = ^ (19) 
Qo 

Above a critical supply of energy, which is expressed in terms of the conversion parameter d 2 , we 
find the occurence of two new solutions for the stationary velocity, corresponding to the active 

(2) (3) 

modes of motion. For F = 0, Vq and Vq both have the same amount but different directions. 
Below the critical threshold the passive mode of motion, = 0, is the only stable solution. At 
the bifurcation point, a change of the stability occurs, consequently above the critical value of d 2 
the active modes of motion are assumed stable. 

For the case F = const. / 0, the analysis of Eq. (17) leads to a different bifurcation diagram 
shown in Fig. lb. Above a critical supply of energy, we find the appearance of two high velocity 
or active modes of motion. One of these active modes has the same direction as the driving force, 
thus it can be understood as the continuation of the normal solution. As Fig. lb shows, the former 
passive normal mode, which holds for subcritical energetic conditions, is transformed into an active 
normal mode, where the particle moves into the same direction, but with a much higher velocity. 
Additionally, in the active mode a new high-velocity motion against the direction of the force F 
becomes possible. While the first active mode would be rather considered as a normal response to the 
force F, the second active mode appears as unnormal (or non-trivial) response, which corresponds 
to an "uphill" motion (cf. Fig. 2). 
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Figure 1: Stationary velocities vq, Eq. (17), vs. conversion rate d^- (a: top) for F = 0, (b: bottom) 
for F = +7/8. Above a critical value of d,2, a negative stationary velocity indicates the posssibility 
to move against the direction of the force. Parameters: qo = 10, 70 = 20, c = 0.01. 

It is obvious that the particle's motion "downhill" is stable, but the same does not necessarily 
apply for the possible solution of an "uphill" motion. Thus, in addition to Eq. (17) which provides 
the values of the stationary solutions, we need a second condition which guarantees the stability of 
these solutions. Before this is carried out, we want to derive a handy expression for the stationary 
velocities in the case F / 0. With the assumption that the term cF is small, the stationary solutions 
of Eq. (17) can be given as: 

^ = 0, v^ = f±Jf s+ (^-^) (20) 
2 7o y 4 7 ^ V70 da) 

We note that Eq. (20) is a sufficient approximation for the stationary velocities, especially in a 
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certain distance from the bifurcation point, thus it will be also used in Sect. 4.1. 

U(x) 




Figure 2: Sketch of the one-dimensional deterministic motion of the particle in the presence of a 
constant force F = — VC/(x) = const.. Provided a supercritical amount of energy from the depot, 
the particle might be able to move "uphill", i.e. against the direction of the force. 



3.2 Stability Analysis 

For the stability analysis, we consider small fluctuations around the stationary values, vq and eQ\ 



v = vq + Sv ; e = eo + Se ; 

Inserting Eq. (21) into Eq. (15), we find after linearization: 

Sv = 
Se = 



Sv 




Se 






eo 



< 1 



7o + d 2 e j + Se \d 2 v 
2d 2 e fo^ + 5e y - c 



With the ansatz: 

Sv ~~ exp {At} ; Se ~ exp {Xt} 
we find from Eq. (22) the following relation for A: 

A (1 ' 2) = --(lo + c + d 2 vl- d 2 e ) 



±y ^(70 + c + d 2 vl - cfoeo) 

c(7o - d 2 e ) - d 2 fo(7o + d 2 e ) 



(21) 



(22) 



(23) 



(24) 



In general, we need to discuss Eq. (24) for the three possible solutions vq which result from Eq. 
(17). Dependent on whether the A for each solution have real or complex positive or negative values, 
we are able to classify the types of the possible stationary solutions in this case. The results are 
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<d 2 < 0.027 : 1 stable node 
0.027 <d 2 < 0.046 : 1 stable node, 1 instable 

node, 1 saddle point 
0.046 < d 2 < 2.720 : 1 stable focal point, 1 

instable focal point, 1 
saddle point 

2.720 < d 2 : 2 stable focal points, 1 

saddle point 



Table 1: Results for the stability analysis, Eq. (24). 

summarized in Table 1. The phase plots shown in Fig. 3a,b present more details. Further, Fig. 4 
shows the real part 9ft(A) of Eq. (24) for the active mode corresponding to the "uphill" motion of 
the particle, which is the most interesting one. 

We find that below the bifurcation point which is = 0.027 for the given set of parameters, only 
one stable node exists in the {v, e} phase space, which corresponds to the passive normal mode. 
Then, a subcritical bifurcation occurs which leads to 3 stationary solutions: a stable and an instable 
node, and a saddle point, since all the A are real. At d 2 = 0.046, however, the nodes turn into focal 
points. With respect to the "uphill motion" we find in Fig. 4 the occurence of an instable node at 
d 2 = 0.027, which then becomes an instable focus for 0.046 < d 2 < 2.720. The respective real parts 
of A are equal in this range, i.e. the A 1 ' 2 are complex. The stability condition is satisfied only if 
R(A) < 0, which is above a second critical value d 2 lt = 2.72 for the given set of parameters. That 
means, for d 2 > 2.72, the instable focal point becomes a stable focus, which is also clearly shown in 
the phase plots of Fig. 3a,b. In both Figures, we see a stable focal point for positive values of the 
velocity, v, which correspond to the stable motion "downhill", i.e. in the direction of the driving 
force. For d 2 = 2.0 < d 2 rit , the phase plot for negative values of v shows an instable focal point, 
which turns into a stable focal point for d 2 = 4.0 > d 2 lt . 

Thus, we can conclude that for d 2 < d 2 an active mode of motion becomes possible, which also 
implies the possibility of an "uphill" motion of the particle. However, only for values d 2 " lt < d 2 , we 
can expect a stable motion against the direction of the force. The interesting result of a possible 
stable "uphill-motion" of particles with internal energy depot will be employed in the following 
section, where we turn to a more sophisticated, piecewise linear potential. 

For our further investigations, it will be useful to have a handy expression for the critical supply of 
energy, d 2 rtt , which allows a stable "uphill" motion. This will be derived in the following, with only 
a few approximations. For the parameters used during the computer simulations discussed later, 
Fig. 4 and Table I indicate that the square root in Eq. (24) is imaginary, thus the stability of the 
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Figure 3: Phase trajectories in the {v, e} phase space for the motion into or against the direction 
of the driving force, which correspond either to positive or negative velocities, (a: top) d,2 = 2.0 
corresponding to an instable "uphill" motion, (b: bottom) 6,2 = 4.0 corresponding to a stable 
"uphill" motion, respectively. Other parameters see Fig. 1. 

solutions depends on the condition: 

7o + c + d2(wg-eo)>0 (25) 
If we insert the stationary value eo, Eq. (16), Eq. (25) leads to a 4th order inequation for vq to 
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Figure 4: Real part of A^ 1 ' 2 ), Eq. (24), vs. conversion parameter d 2 for the stationary motion against 
the force, F. Parameters see Fig. 1. 



obtain stability: 



(70c - d 2 q ) < Vq d\ + Vq (j d 2 + 2cd 2 ) + c z 



(26) 



For a stable stationary motion of the particle, both Eq. (17) and Eq. (26) have to be satisfied. 

The critical condition for stability just results from the equality in Eq. (26), which then provides 
a replacement for the prefactor (70c — d 2 qo) in Eq. (17). If we insert the critical condition into Eq. 
(17), we arrive at a 5th order equation for vq: 



cF 



(27) 



In order to simplify the further discussion, we assume that the internal dissipation is negligible, 
c = 0. Then, Eq. (27) gives the simple nontrivial solution: 



■— ; if c = 
d 2 



(28) 



This expression can be used to eliminate the stationary velocity, vq, in Eq. (26). With the assump- 
tion c = 0, we obtain now from the critical condition, i.e. from the equality in Eq. (26), a relation 
between the force, F, and the conversion parameter, d 2 . Combining Eq. (26) and Eq. (28) results 
in 



(~F)^d 2 + 7o(-^) 2/3 df - q 4 /3 = 



(29) 
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Because of c?2 > 0, the trivial and the negative solution of Eq. (29) can be neglected, and we finally 
arrive at the following critical relation for d,2(F): 

<3o> 

Fig. 5 shows c#f 4i as a function of the force, F. In the limit of negligible internal dissipation, the 
relation d^ lt {F) describes how much power has to be supplied by the internal energy depot in order 
to allow a stable motion of the particle in both directions, in particular a stable uphill motion of 
the particle. 

12 
10 
8 

CM £ 

4 

2 


0.0 1.0 2.0 3.0 

|F| 

Figure 5: Critical conversion rate, d!f 1 *, Eq. (30), vs. amount of the driving force, to allow a 
stable motion of the particle both "downhill" and "uphill" (cf. Fig. 2). For c?2 > the particle 
is able to move also against the direction of the force. Parameters see Fig. 1. 
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4 Deterministic Motion in a Ratchet Potential 



4.1 Simulation Results for a Single Pumped Particle 

For further investigations of the motion of pumped particles, we specify the potential U(x) as a 
piecewise linear, asymmetric potential (cf. Fig. 6), which is known as a ratchet potential: 



U{x) 



^{x - nL} 

if nL < x < nL + b 
j^{(n + l)L-x} 

ifnL + b<x<(n + l)L 
(n = 0,1,2, ... 



(31) 




Figure 6: Sketch of the asymmetric potential U(x) (eq. (31)). For the computer simulations, the 
following values are used: b=4, L=12, Uq = 7 in arbitrary units. 

Further, we will use the following abbreviations with respect to the potential U(x), Eq. (31). The 
index i = {1,2} refers to the two pieces of the potential, l± = b, h = L — b. The asymmetry 
parameter a should describe the ratio of the two pieces, and F = —VU = const, is the force 
resulting from the gradient of the piecewise linear potential. Hence, for the potential U(x), Eq. 
(31), the following relations yield: 



Fi = 
Fi = 



^; Fo 



b 
Uq 
L 



U 



a = — 



(1 + a) ; F 2 = 



h 

b ' " h 
UqI + o 



L-b 



>2 



L 



(32) 



Whether or not the particle will be able to leave one of the potential wells described by Eq. (31), 
depends in a first approximation on the height of the potential barrier, Uq, and on the kinetic 
energy of the particle. For particles with internal energy depot, the actual velocity depends also 
on the conversion of internal into kinetic energy, Eq. (7). In agreement with the investigations in 
the previous section, we can in principle distinguish between two different types of motion: (i) a 
bound motion, i.e. the particle will not leave the potential well because of the subcritical supply 
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of energy from the depot, but its position might oscillate within the boundaries, (ii) an unbound 
motion, i.e. the particle will be able to leave the potential well because of the supercritical supply 
of energy and move freely. Both types of motion have analogies to the localized and delocalized 
states of electrons in solid state physics as will be discussed in more detail in a forthcoming paper. 

In the following, we discuss computer simulations of the deterministic motion of one pumped 
Brownian particle in a ratchet potential. The particle (mass m = 1) starts its motion outside the 
potential extrema; hence, there is an initial force on the particle. The results for a single particle 
are shown in Fig. 7, where two different sets of parameters are used: 

(i) a small internal dissipation, c, which means a nearly ideal energy depot, and a large friction 
coefficient, 70, resulting in a strongly damped motion, 

(ii) an internal dissipation, c, 10 times larger, an energy influx, qo, ten times smaller, and a friction 
coefficient, 70, 100 times smaller than in (i), resulting in a weakly damped motion, 

We note, that in the computer simulations always the complete set of Eq. (7) for the particles 
is solved, regardless of the possible approximations. The trajectories x(t) in Fig. 7 indicate a 
nearly uniform motion of the particle into one direction. The continuous motion corresponds to a 
delocalized state of the particle in the ratchet potential, whereas for subcritical energetic conditions 
only localized states exist. The transition from localized to delocalized states will be discussed in 
more detail in a forthcoming paper [31]. 

As shown in Fig. 7a, the less damped motion of a single particle may result in steady oscillations in 
the velocity, v(t), and the energy depot, e{t). Only if the damping is large enough, the velocity and 
the energy depot may reach constant values (cf. Fig. 7b). These values are of course different for 
each piece of the potential, hence the periodical movement through every maximum or minimum 
of the potential results in jumps both in the velocity and the energy depot, which are followed by 
oscillations. In the phase space shown in Fig. 8a, the motion of the particle appears as a transition 
between two stable fix points, each of which describes the stable motion on one flank. 

In the following, we restrict the discussion to the strongly damped case. The oscillations which occur 
in v and e are damped out on a characteristic time scale, r = 1/70- If we assume that the particle 
moves on the different pieces of the potential {b, L — b} (cf. Fig. 6) during the two characteristic 
time intervals: T& = b/vb and = (L — then the particle is subject to a constant force 

only as long as r C or r < Th-h, respectively. For times larger than the characteristic time, r, 
the motion of the particle can be described by the equation of the overdamped limit, Eq. (11). If 
we neglect again stochastic influences, Eq. (11) can be rewritten in the form: 

= -[70-cfeen] v + F (33) 

where |a;| = \vq\ = const, is the velocity in the overdamped limit and F = {iq,^} is defined by 
Eq. (32). The stationary value for the internal energy depot, eo, is given by Eq. (16). If we assume 
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Figure 7: Trajectory x(t), velocity v(t) and energy depot e(t) for a single particle moving in a 
ratchet potential (Fig. 6). Parameters: (a: top) qo = 1.0,7 = 0.2, c = 0.1, ofo = 14.0, (b: bottom) 
go = 10, 7o = 20, c = 0.01, d 2 = 10. Initial conditions: x(0) G [4, 12], u(0) = 0, e(0) = 
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Figure 8: (a: top) Energy depot e{t) vs. velocity v(t), (b: bottom) prefactor: 7 — d,2e(t) for the 
strongly damped motion of a single particle (data from Fig. 7b). 

again that the term cF is small, then the constant velocity can be calculated from Eq. (20), with 
F specified as F\ or F2, respectively. For any constant force F, there are two possible non-trivial 
solutions of Eq. (20): a positive and a negative velocity with different, but constant amount, which 
depend on the gradient of the potential. The nontrivial values v 7^ of Eq. (20) can be compared 
with the constant values obtained in the simulations, and we find: 

Fig. 7b Eq. (20) 

lower value v 0.665 0.6642 (34) 

upper value v 0.728 0.7288 

The eqs. (11), (33) for the overdamped limit indicate that the dynamics remarkably depends on the 
sign of the prefactor 70 — cfeeo, which governs the influence of the potential and the stochastic force. 
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Therefore, the prefactor should be discussed in more detail, now. Fig. 8b shows that the prefactor 
7o — (faeit) displays a behavior similar to the velocity, Fig. 7b. The prefactor jumps between a 
positive and a negative constant value, which can be approximated by means of the constant, en, 
Eq. (16), reached after the oscillations damped out. It is shown that the jump occurs at the same 
time when the gradient of the potential changes its sign. This can be also proved analytically. Using 
Eqs. (10), (20), the prefactor 70 — cfoeo in Eq. (11) and Eq. (33), respectively, can be rewritten and 
we find after a short calculation: 



1 1 



70 - d,2 eo 270 Fi 



Fi± x Ff+A{q 0l0 - C1 l/d: 



(35) 



This means that the product of the prefactor and the potential gradient always has the same 
positive (or negative) sign, and the direction of motion for the particle is only determined by the 
initial condition. 

The prefactor 70 — c^en describes the balance between dissipation and the energy supply from the 
internal depot of the particle. Therefore it is expected that the time average of the prefactor should 
be zero, when averaged over one time period T: 

(70 - d 2 e ) T = (36) 

This can be proven using T = T5 + as discussed above. The values vj,, vi-b and e&, e^-b for 
the related velocity and energy depot along the different pieces of the potential can be calculated 
by means of Eqs. (20), (33) which leads directly to Eq. (36). 



4.2 Investigation of the Net Current for an Ensemble of Particles 

Let us now discuss the deterministic motion of an ensemble of N pumped Brownian particles in a 
ratchet potential. For the computer simulations, we have assumed that the start locations of the 
particles are equally distributed over the first period of the potential, {0, L} and their initial velocity 
is zero. In the deterministic case, the direction of motion and the velocity at any time t are mainly 
determined by the initial conditions. Hence, particles with an initial position between {0, b}, which 
initially feel a force into the negative direction, most likely move with a negative velocity, whereas 
particles with an initial position between {b, L} most likely move into the positive direction. This 
is also shown in Fig. 9a, where the velocity v is plotted versus the initial position of the particles. 
Oscillations occur only at the minima and maxima of the related potential, indicating a strong 
sensitivity to the initial condition in these regions. The distribution of the final velocity is shown 
in Fig. 9b. 

From Fig. 9b we see two main currents of particles occuring, one with a positive and one with a 
negative velocity, which can be approximated by Eq. (20). The net current, however, has a positive 
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Figure 9: (a: top) Final velocity v e after t=2.000 simulation steps (averaged for 10.000 particles) vs. 
initial location xq of the particles, (b: bottom) Distribution of the final velocity v e . Parameters see 
Fig. 7b, initial locations of the particles are equally distributed over the first period of the ratchet 
potential {0, L}. 

direction, since most of the particles start with the matching initial condition. The time dependence 
of the averages is shown in Fig. 10. The long-term oscillations in the average velocity and the average 
energy depot result from the superposition of the velocities, which are sharply peaked around the 
two dominating values (cf. Fig. 9b). 

The existence of periodic stationary solutions, vq(x) = vq(x ± L), requires that the particles are 
able to escape from the initial period of the potential. For a continuous motion which corresponds to 
delocalized states, the particles must be able to move "uphill" on one or both flanks of the ratchet 
potential. In Sect. 3, we already investigated the necessary conditions for such a motion for a single 
flank, and found a critical condition for the conversion rate, cfo, Eq. (30). In order to demonstrate 



19/27 



Frank Schweitzer, Benno Tilch, Werner Ebeling: 
Uphill Motion of Active Brownian Particles in Piecewise Linear Potentials 
European Journal of Physics B (1999, in press) 




Figure 10: Averaged location (x), velocity (v) and energy depot (e) of 10.000 particles vs. time t. 
Parameters see Fig. 7b. 

the applicability of Eq. (30) for the ratchet potential, we have investigated the dependence of the 
net current, expressed by the mean velocity (v), on the conversion rate, cfe, for the overdamped 
case. The results of computer simulations are shown in Fig. 11. 

In Fig. 11, we see the existence of two different critical values for the parameter d%, which correspond 
to the onset of a negative net current at d^"* 1 an d a positive net current at d?f tt2 . For values of cfe 
near zero and less than d" 111 , there is no net current at all. This is due to the subcritical supply 
of energy from the internal depot, which does not allow an uphill motion on any flank of the 
potential. Consequently, after the initial downhill motion, all particles come to rest in the minima 
of the ratchet potential, with vq = as the only stationary solution for the velocity. With an 
increasing value of d2, we see the occurence of a negative net current at dSf** 1 . That means, the 
energy depot provides enough energy for the uphill motion along the flank with the lower slope, 
which, in our example, is the one with F = 7/8 (cf. Fig. 6). If we insert this value for F into the 
critical condition, Eq. (30), a value df 111 = 2.715 is obtained, which agrees with the onset of the 
negative current in the computer simulations, Fig. 11. 
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Figure 11: Average velocity (v) vs. conversion parameter di- The data points are obtained from 
simulations of 10.000 particles with arbitrary initial positions in the first period of the ratchet 
potential. Other parameters see Fig. 7b. 

For drj™* 1 < d,2 < d?, 7 " 112 , a stable motion of the particles up and down the flank with the lower slope is 
possible, but the same does not necessarily apply for the steeper slope. Hence, particles which start 
on the lower slope with a positive velocity, cannot continue their motion into the positive direction 
since they are not able to climb up the steeper slope. Consequently, they turn their direction on 
the steeper slope, then move downhill driven by the force into the negative direction, and continue 
to move into the negative direction while climbing up the lower slope. Therefore, for values of the 
conversion rate between dg™* 1 and c^"* 2 , we only have an unimodal distribution of the velocity 
centered around the negative value: 



which is independent of d2 if the term (c/d2) in Eq. (20) is negligible, which holds for the considered 
case. For d2 > d?f it2 , the energy depot also supplies enough energy for the particles to climb up the 
steeper slope, consequently a periodic motion of the particles into the positive direction becomes 
possible, now. In our example, the steeper slope corresponds to the force F = —7/4 (cf. Fig. 6) 
which yields a critical value d^ ltl = 5.985, obtained by means of Eq. (30). This result agrees with 
the onset of the positive current in the computer simulations, Fig. 11. 

For c?2 > c#f rf2 , we have a bimodal velocity distribution, as also shown in Fig. 9b. The net current, 
which results from the average of the two main currents, has a positive direction in the deterministic 
case, because most of the particles start into a positive direction, as discussed above. We may simply 
assume, that the number of particles in each direction is roughly proportional to the length of the 




(37) 
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flank where they started from, which is also indicated by the velocity distribution, Fig. 9b. Then 
the mean velocity in the strongly damped case can be approximated by: 



1 N 1 



3 «i + 3^ 



(38) 



where v\ and v 2 are the stationary velocities on each flank, which, in the limit of an nearly ideal 
energy depot, can be determined from Eq. (20). With the negative velocity, v\, Eq. (37), and the 
positive velocity, 



V2 



+ 



F 2 



270 V 4 ^o 



+ 



go 

To 



0.664 for F 



7 

4 



(39) 



we find from Eq. (38) for d 2 > d^ %t2 an average velocity, (v) = 0.216, which also agrees with the 
computer simulations, Fig. 11. 

The results of the computer simulations have demonstrated that in the deterministic case for the 
given special initial condition, i.e. the equal distribution of particles over the ratchet period, the 
direction of the net current can be adjusted by choosing the appropriate values of the conversion 
rate, d 2 . 

The critical values for d 2 , on the other hand, depend on the slope of the two flanks of the potential, 
expressed by the force F. Lower slopes also correspond to lower values of the conversion rate, 
because less power is needed for the uphill motion. 

We conclude our results by investigating the influence of the slope on the establishment of a positive 
or negative net current. With a fixed height of the potential barrier, Uq, and a fixed length L, the 
ratio of the two different slopes is described by the asymmetry parameter a = l 2 /l\ = —F1/F2, Eq. 
(32). The occurence of a current in the ratchet potential requires the possibility of uphill-motion, 
which depends on the critical supply of energy, described by Eq. (30). In order to obtain the critical 
value for the asymmetry of the potential, we replace the force F in Eq. (30) by the parameter a, 
Eq. (32). In our example, the flank l\ of the potential has the steeper slope, so the critical condition 
is determined by F± = Uq/L (1 + a). As the result, we find: 



crit 



L 



2 2 + 



— d 2 +q d 2 ' 



3/2 



(40) 



a > 1 gives the critical value for the asymmetry, which may result in a reversal of the net current. 
For a > a crtt , the flank l\ is too steep for the particles, therefore only a negative current can occur 
which corresponds to the unimodal velocity distribution discussed above. For 1 < a < a crtt , however, 
the particles are able to move uphill either flank. Hence, also a positive current can establish and 
the velocity distribution becomes bimodal, which results in a positive net current. 
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Figure 12: Average velocity (v) vs. asymmetry parameter a, Eq. (32). The data points are obtained 
from simulations of 10.000 particles with arbitrary initial positions in the first period of the ratchet 
potential, (o): d 2 = 10 ( ), (O): d.2 = 20 ( ). Other parameters see Fig. 7b. 

The current reversal from a negative to a positive net current is shown in Fig. 12. Dependent on the 
value of the conversion rate, cfe, we see the switch from the negative to the positive value of the net 
current at a critical value of the asymmetry parameter a. Because of the definition of a, the results 
for a < 1, are the inverse of the results for a > 1. Obviously, for a symmetric ratchet, a = 1 no net 
current occurs, because the two main currents compensate. From Eq. (40), we obtain a crtt = 3.5 for 
c?2 = 10, and a cnt = 6.5 for d,2 = 20, which both agree with the results of the computer simulations, 
Fig. 12. Further, the results of Fig. 12 shows that the stationary velocities are independent of d-2 
in the limit of an nearly ideal energy depot, which is also indicated by Eq. (20). 



5 Summary and Conclusions 

In this paper, we have investigated the dynamics of Brownian particles with an internal energy 
depot, which can be filled by the take-up of energy from the environment. This extension which is 
inspired by the biological features of active motion, has several aspects: (i) the spatial or spatial- 
temporal inhomogeneous distribution of energy from the environment can be considered, (ii) the 
supply of depot energy for different activities, i.e. for accelerated motion or signal-response behav- 
ior, can be modeled, (iii) in addition to external dissipative processes caused by friction, internal 
dissipation can be included, (iv) the adjustment of the time scale for the relaxation of the energy 
depot allows to consider delay effects, i.e. for the acceleration of motion. 

In this paper, we are mainly interested in the question how the internal energy depot will affect 
the particle's motion in an external potential. We found that provided some supercritical energetic 



23/27 



Frank Schweitzer, Benno Tilch, Werner Ebeling: 
Uphill Motion of Active Brownian Particles in Piecewise Linear Potentials 
European Journal of Physics B (1999, in press) 



conditions, the Brownian particles are able to move in a "high velocity" or active mode of motion 
characterized by a velocity much larger than the Stokes velocity. The latter one can be denoted 
as the passive mode of motion. In addition to stochastic forces, it is basically determined by the 
response to the gradient of an external potential. One of the possible active modes of motion can be 
understood as a continuation of the passive or Stokes mode, the particles motion in the direction 
of the force being accelerated. 

Additionally, we also found an active mode which describes a motion of the particle against the 
gradient of an external potential. We have investigated the critical conditions for the deterministic 
motion of the particle in the overdamped limit for the case of a linear potential. We found that 
the "uphill" motion is described by two critical conversion rates. The bifurcation diagram shows 
that, at a critical value d^ , the possibility of an "uphill" motion appears as a new solution for the 
stationary velocity, i.e. "uphill" motion requires the existence of multiple stationary solution for 
the velocity. However, this active mode remains unstable as long as the conversion rate is below a 
second critical value dZf 1 *, Eq. (30), which depends on parameters describing the energy balance, 
such as the friction coefficient, 70, and the influx of energy into the depot, qo, and on the gradient 
of the potential, i.e. the resulting force F. A stable "uphill" motion occurs only for a supercritical 
conversion rate. 

These results allow us to interpret the deterministic motion of an ensemble of particles in a ratchet 
potential, i.e. a piecewise linear, periodic asymmetric potential with two different slopes. Initially, 
the particles are equally distributed over the first period of the potential. In order to produce a 
current, the particles must be able to escape from the potential well, which depends on the supply 
of energy. We are able to find two different critical conversion rates, which describe the onset of a 
directed current into the different directions. As long as the particles are only able to move uphill the 
flank with the lower slope, we have an unimodal velocity distribution, corresponding to a negative 
net current. But if the particles are able to move uphill either flank, we find a bimodal velocity 
distribution, and a positive net current for the given conditions. Hence, in the deterministic case 
the direction of the net current can be controlled via a single parameter, di, which describes the 
conversion of internal into kinetic energy. We are able to calculate the critical parameter, as well 
as the resulting velocity of the net current from analytical considerations, which agree with the 
results found in the computer simulations. Alternatively, the critical supply of energy can be also 
reformulated in terms of a critical asymmetry of the ratchet potential, which restricts the occurence 
of a (positive or negative) net current. 

To conclude the results obtained, we have shown that an ensemble of pumped Brownian particles 
moving in a ratchet potential can produce a directed net current. In this respect, our result agrees 
with the conclusions of other physical ratchet models which have been proposed to reveal the 
microscopic mechanisms resulting in directed movement. Due to [19, p. 295] a ratchet system which 
is meant to be "a system that is able to transport particles in a periodic structure with nonzero 
macroscopic velocity although on average no macroscopic force is acting." Indeed, we have shown 
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in Eq. (36) that for the stationary approximation (70 — cfeeo) t = holds, i.e. the force acting 
on the particle is of zero average with respect to one period, r. However, different from other 
ratchet models, ie. the rocking ratchet [27, 28, 30] or the diffusion ratchet [29] which assume a 
spatially uniform time-periodic force, the force in our model switches between two constant values, 
dependent on the moving direction and the flank the particle is moving on. Hence, it does not 
represent a spatially uniform force, and we may conclude that the mechanism of motion which Eq. 
(11) is based on, should be different from the previous mechanisms which originate directed motion 
in a ratchet potential. 

In this paper, we have restricted the discussion to the deterministic motion of Brownian particles 
with an internal energy depot. We note that the stochastic motion of an ensemble of Brownian 
particles in a ratchet potential will be investigated in detail in a forthcoming paper [31], which pays 
particular attention to the influence of stochastic effects on the establishment of the net current. 
Additionally, we will also discuss the non-equilibrium distribution functions for Brownian particles 
with internal energy depot [26]. 

As a final remark, we note in a more general sense that the introduction of an internal energy 
depot adds an interesting and new element to the known model of Brownian particles, which can 
be useful in two different respects. With respect to the physical aspects, the dynamical system now 
has a new degree of freedom, which increases the phase space, V = {xi,Vi, e}. Due to the additional 
energy supply, the system is driven into non-equilibrium. Hence, a qualitative new behavior in the 
particle's motion can be obtained. 

With respect to possible biological aspects, the extension of the Brownian particle model by mech- 
anisms of energy take-up, storage and conversion, should contribute to the development of a mi- 
croscopic theory of active biological motion. The final goal of such a project could be a microscopic 
image of well known phenomenological models of biological motion, taking into account energy 
balances that are related to the mechanisms of energy pumping and energy dissipation. 
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